Possible Evidence for A Slush Vortex Phase in the Vicinity of Vortex Glass Transition 

in YBa 2 Cu 3 7 5 Thin Films 
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Current-voltage I — V curves have been carefully measured for YBa2Cu307_a thin films by following 
different thermal or field annealing procedures. Although all data can be described quite well by 
the vortex-glass theory, it is found that the dissipation in a small region just above the vortex- 
glass transition is history dependent: the dissipation in the warming up or field-increasing process is 
larger than that in the cooling down or field-decreasing process. We attribute this history dependent 
dissipation to possible existence of a slush vortex phase. 
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The vortex matter continues to be an attractive sub- 
ject in the study of high temperature superconductors 
(HTS). Due to the intrinsic properties {e.g., small co- 
herence length £, high operating temperature T, high 
anisotropy e, low charge density or high normal state re- 
sistivity p n , etc. ) the HTSs have a rich and interesting 
vortex phase diagram which has received enormous ef- 
fort in the past decade Much of the complexity of 
the vortex dynamics arises from the competing roles of 
vortex-vortex interactions, random pinning, and thermal 
fluctuations. In the absence of quenched disorder, the 
ideal vortex matter should have a crystalline ground state 
and can undergo a genuine solid-liquid phase transition 
upon heating. In real type-II superconductors, however, 
the long-range translational order of the vortex lattice 
is destroyed by random impurities Fisher, Fisher, 
and Huse || conjectured that the low-temperature vor- 
tex phase in systems with random pinning could be a vor- 
tex glass (VG) characterized by the extremely non-linear 
V(I) (or E(j)) dissipation p|, or E ~ exp[— (j c /j)1, and 
non-logarithmic magnetic relaxatio n HL consistent with 
many of the experimental findings n-S]. In HTS films 
with strong random pinning centers, early experiments 
have shown a continuous transition with no hysteresis in 
transport properties ||. In this Letter, we report new 
results on similar HTS thin film samples using a more 
careful measurement procedure. Our data show clear ev- 
idence of a history dependent dissipation in the vicinity 
of the VG transition. 

Thin films measured for this work were prepared by 
magnetron sputtering deposition on (100) MgO sub- 
strates. The crystallinity of the films has been checked 
by the X-ray diffraction (XRD) pattern which shows only 
the (001) peaks indicating a good epitaxy. The oxygen 
in the films has been depleted by annealing them in low 
pressure oxygen atmosphere, therefore the samples are 
slightly underdoped. The films with thickness of about 
2000 A have been patterned into a bridge shape with 



width of 4 /im and length of 10 fim. After patterning, 
four silver pads were deposited and a post annealing at 
about 400 °C was made to reduce the contact resistance. 
(Joule heating is negligible for the current range of inter- 
est here.) The final bridge has a zero resistance temper- 
ature T c o = 86K and transition width less than 1 K, the 
critical current density j c {77K, H = 0) = 2 x 10 6 A/cm 2 
using a voltage criterion of E c = 1 /iV/mm. The trans- 
port measurement has been carried out with a Keithley 
220 programable current source and a 182 nano-volt mul- 
timeter. An Oxford cryogenic system with magnetic field 
up to 8 Tesla was used to vary the temperature and mag- 
netic field. The temperature is controlled by adjusting 
the PID parameters of Oxford ITC4 tempearture con- 
troller with a temperature accuracy of about 0.02 K. 

Fig. 1 shows the resistive transitions under several 
magnetic fields. It is found that the R(T) curves shift 
parallel to lower temperatures being similar to that in 
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FIG. 1. Resistive transitions at magnetic fields ranging 
from T to 8 T. The parallel shift of R(T) is a property due 
to underdoping. Inset shows the upper critical field H C 2(T) 
and the vortex glass transition line H g (T). 
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the underdoped region as observed by Suzuki et al. |T(]] 
and Wen et al. Q This parallel shift of R(T) curve 
may be connected with the electronic properties in un- 
derdoped region, therefore we leave the discussion to a 
separate publication. The upper-critical field H C 2 (T) can 
be determined by extrapolating the flat line in the normal 
state and the steep transition portion to a cross point. 
As shown in the inset to Fig.l, the H C 2(T) curve has 
a slope of about 0.8 T/K near T c (0). The irreversibil- 
ity transition temperature Ti rr was determined by taking 
the criterion of R = 1CP 4 J7. It is found that Ti rr (H) is 
close to the VG transition temperatures T g determined 
by extrapolating the curve of dT/dlnR vs. T to zero ( 
dT/dlnR = ). Between H c2 (T) and H g (T) it is sup- 
posed to be a 3D vortex liquid state separated from the 
VG state by the H g (T) line. 

I — V curves have been carefully measured at mag- 
netic fields of 0.5 T, 1 T, 2 T, 4 T and 8 T. For mea- 
suring these curves, four different annealing procedures 
are used: warming up (T up ) and cooling down (Tdown) 
with fixed fields; field ascending (H up ) and field descend- 
ing (Hdown) & t fixed temperatures. Special cautions were 
taken in these processes. In the T up (Tdown) process, the 
sample is first cooled in a magnetic field to a temperature 
which is 10 K below (above) the target temperature, then 
it is heated up (cooled down) by stablizing the temper- 
atures in steps of 1 K up (down) to 1 K below (above) 
the target value. The next step is to approach the target 
temperature with an extremely slow speed by stablizing 
T in steps of 0.1 K. In this process, over-heating or over- 
cooling should be avoided in preparing a reliable vortex 
state which "records" the history. In the H up (Hdown) 
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FIG. 2. I — V curves measured in the so-called H up and 
Hdown processes (see text) from 76 K to 88 K with increments 
of 0.3 K. It is found that, in the temperature region just above 
T g ( marked here by the thick dotted line ) the dissipation 
in the H up process is always higher than that in the Hdown 
process showing a history dependent effect. It disappears in 
the high current (voltage) region for all temperatures. 



process, the temperature is first fixed, the field is swept 
up (down) to the target field with a sweep rate of about 
100 G/s. The same cautions should be taken in order to 
prevent the field over-shooting. 

Fig. 2 presents the I — V curves measured in H up (solid 
line) and Hdown (dashed line) processes. It is found that 
at temperatures below and well above T g marked here by 
the thick dotted line, the two sets of data coincide very 
well. In the small region just above T g , it is surprising to 
see that, the I— V curves measured in the H up process lies 
always above that measured in Hdown process indicating 
a history dependent dissipaption. When the current is 
increased into the region with a common negative curva- 
ture of InV vs. Inj, this effect disappears gradually. For 
a closer inspection to this history dependent dissipation, 
an enlarged view in the voltage region between 5 x 10~ 6 
V and 1 x 10 -4 V is shown in Fig. 3. This effect has been 
observed on the same sample with different runs of mea- 
surement by either varying the magnetic field at a fixed 
temperature or varying temperatueres at a fixed mag- 
netic field. It is also repeatable on other YBa2Cu307_5 
thin films. Therefore it is a real effect. 

We find that, in spite of the striking history effects, 
the I — V curves can nevertheless be shown to satisfy the 
scaling relations predicted by the VG theory. According 
to Fisher et al. [3] , at a second-order phase transition the 
I — V curves can be scaled according to 



j| 7 1 _ j 1 ^(z+2-d) T \ T — T 

where v and z are critical exponents characterizing the 
VG transition, d is the dimensionality and f(x) is a uni- 
versal function for all temperatures at a given magnetic 
field. In Fig. 4 we present the VG scaling according to 
eq.(l) for the I — V curves measured in the H up process 
at 1 T. Since the critical region for a VG scaling is 
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FIG. 3. An enlarged view of Fig.2 in the voltage region 
between 5 x 10~ 6 V and 1 x 10~ 4 V. The history dependent 
effect is more clear here. 
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narrow, therefore we use the data from 79 K to 85 K 
and voltage below 5 x 10~ 3 V for the scaling. As shown 
by Fig. 4, the data can indeed be scaled onto two main 
branches. The glass exponents obtained here are v = 1.5 
and z = 5 which are very close to the values found by 
Koch et al. Q in the same type of films. Another check is 
the ohmic resistivity just above T g which, in VG model, 
reads p — po{T — T g ) s with s = v(z — d + 2), where d 
is the dimensionality taken to be 3 for our case. Taking 
v = 1.5 and z = 5 obtained from the scaling it is found 
that s = 6. In the inset to Fig. 4 we show the data of 
ohmic R(T). The s value derived from the R(T) data 
is 6.1 which is almost identical to that obtained from 
the scaling. The VG scalings for the field of 2 T, 4 T 
and 8 T give almost the same glass exponents showing 
the generality of our films. (It is important to point out 
that the apparent criticality in the I — V characteristics 
could be caused by the correlated disorder such as grain 
boundaries and twins. In that case, the universality class 
should be that of a Bose glass |pq| . For the purpose of 
comparison, we followed the same procedures as that of 
Koch et al. §) 

To further investigate the history dependent effect, we 
take the difference between the currents corresponding 
to the same voltage in the increasing field (temperature) 
and deceasing field (temperature) processes at fields of 1 
T and 4 T. As shown in Fig. 5 (a) to (d), a clear peak 
appears in a narrow region just above the VG transition 
temperature T g marked by vertical thick dashed lines for 
all measurements. The open symbols: squares, circles 
and up-triangles represent the data taking at voltages of 
5 x 10" 6 V, 5 x 1(T 5 V and 5 x 1(T 4 V, respectively. 
Clearly, the peak height due to this history dependent 
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FIG. 4. Vortex glass scaling of the data shown in Fig. 2. 
The same scaling has been taken to the data measured at 
other fields 1 T, 2 T, 4 T, and 8 T all yielding the similar 
values for v and z. Inset shows the relation of resistance 
versus T — T g from which it is found that s — 6.1 being very 
close to that determined from the scaling. 



dissipation is supressed at a high voltage or current 
where the logl vs. logV shows a common negative cur- 
vature. This indicates that the history dependent ef- 
fect appears in the so-called thermally-activated-flux-flow 
(TAFF) regime (ohmic). Another interesting finding is 
that the hysteretic behavior is stronger in the process of 
varying temperature than varying magnetic field. 

The history dependent effect in transport properties 
has been observed in weak-pinning YBa2Cu307~,5 near 
the first-order phase transition O], In thin films with 
much stronger disorder, the transition is supposed to 
be second order |3| and any hysteretic behavior is unex- 
pected. Our data show that however even in thin films, 
the hysteretic behavior, although weak, can still be ob- 
served. This hysteresis may be understood by assuming 
the existence of a slush vortex phase just above T g . The 
main argument is that the vortex matter with this hys- 
teretic behavior is not in a single liquid state, rather it is 
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FIG. 5. The difference between the currents correspond- 
ing to the same voltage in the increasing field (temperature) 
and deceasing field (temperature) processes at fields of 1 T 
and 4 T. A clear peak appears in the small region just above 
T g marked here by the vertical thick dashed line. The open 
symbols: squares, circles and up-triangles represent the data 
taking at voltages of 5 x 10" 6 V, 5 x 1(T 5 V and 5 x 10" 4 V, 
respectively. In the large voltage region, the peak is strongly 
suppressed. 
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When T is below T g , the vortex system is in the glassy 
state but with finite shear rigidity, therefore the small 
droplet or crystallites with translational order will play 
a role as building blocks for the whole elastic VG state. 
This picture is also consistent with the collective creep 
theory Q that at a current smaller than the critical value, 
the optimal hopping object will be constructed by many 
Larkin domains. This will lead to a divergent high ac- 
tivation barrier in the small current limit. Above T g , 
the strong thermal fluctuations will destroy the long- 
range phase coherence, leading to a second order melting. 
While just above the glass melting temperature, some 
small droplets or crystallites with short range order may 
survive. The melting of these small vortex cystallites may 
be first order and occur at a temperature slighly higher 
than T g causing the hysteresis effect. Since the vortices in 
the liquid state are softer than those in the small vortex 
crystals, they can easily adapt their configuration to find 
some strong pinning sites. Therefore when the tempera- 
ture is lowered to the small region above T g the vortices 
are still in the plastic state, while when the tempera- 
ture is increased to the same region from below T g some 
small vortex crystals remain unmelted. This will give rise 
to the slightly higher dissipation in the increasing tem- 
perature (field) process than in the decreasing process. 
Computer simulations by Dominguez et al. @ have in- 
dicated that a large current can promote the melting of 
vortex lattice. This may explain in our present case why 
in the high current part the history dependent effect dis- 
appears. 

The term "slush vortex" was first introduced by Wor- 
thington et al. |l5| ] in explaining the two step resistive 
transition in the defect-enhanced YBa2Cu307_^ single 
crystals. Recently Nishizaki et al. |l6| have found also in 
underdoped YBa2Cu307_a single crystals that the vor- 
tex glass transition T g (H) can be well below the crit- 
ical point of melting line T m (H) showing the possible 
existence of a vortex slush regime between T g (H) and 
T m (H). Above two cases are close to the clean limit with 
single crystals as the beginning samples. Our present 
work reveals however that this slush vortex phase may 
also exist in exetremely disordered system, such as thin 
films. It would be interesting to know whether the his- 
tory dependent dissipation found in our present paper 
occurs also in other systems, such as E^S^CaC^Os or 
Tl2Ba2CaCu208 since it has been shown that for 
such highly anisotropic systems, when the magnetic field 
is high enough, the vortex system melts at zero K with 
a non-zero linear resistivity at any finite temperatures. 

In conclusion, by following different thermal and field 
annealing processes, a history dependent dissipation is 
observed in the vicinity of the vortex glass melting. This 
effect is attributed to the possible existence of a slush 
vortex phase due to the non-uniform distribution of the 
pinning sites. 
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